Metal-insulator transition in a 2D electron gas: Equivalence of two approaches for 

determining the critical point 
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The critical electron density for the metal-insulator transition in a two-dimensional electron gas 
can be determined by two distinct methods: (i) a sign change of the temperature derivative of 
the resistance, and (ii) vanishing activation energy and vanishing nonlinearity of current-voltage 
characteristics as extrapolated from the insulating side. We find that in zero magnetic field (but not 
in the presence of a parallel magnetic field), both methods give equivalent results, adding support 
to the existence of a true zero-field metal-insulator transition. 
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The observation of a metal-insulator transition (MIT) 
in two-dimensional (2D) electron and hole systems (see a 
review [n] and references therein) challenges the scaling 
theory of localization which predicts all electron states 
to be localized in an infinite disordered 2D system at 
zero temperature and in zero magnetic field 0. In fact, 
the predicted localization of all carriers is beyond ex- 
perimental verification since in real 2D systems, it can 
break down because of any perturbation such as a finite 
temperature, finite system size, magnetic impurities, etc. 
To estimate the limiting behavior at zero temperature 
experimentally, one can follow two principal approaches: 
extrapolate the temperature dependence of the resistance 
to T = 0, or analyze temperature-independent character- 
istics {e.g., the correlation length). Moreover, compari- 
son of distinct criteria is the experimental test to find out 
whether or not a true MIT exists. The resistance of high- 
mobility silicon and p-GaAs samples is "flat" at a certain 
carrier density in a wide range of temperatures P,H,p|, 
and it is therefore tempting to extrapolate it to T = 
and identify the corresponding electron density with the 
critical density for the metal-insulator transition. How- 
ever, in some of the suggested explanations for the MIT 
(see, e.g., Ref. Iq|), it was stated that within their mod- 
els, the "critical" curve is not flat. Therefore, to verify 
whether or not the temperature-independent resistance 
corresponds to the critical density, an independent de- 
termination of the critical point is necessary. This is why 
in this paper, we compare the critical density obtained 
by a finite-temperature criterion with that obtained by 
a temperature-independent one. Equivalent results from 
these two distinct methods would justify the extrapola- 
tion of the flat resistivity curve to T = and strongly 
support the existence of a true MIT in zero magnetic 
field. 

Methods for determining the MIT point have been de- 
scribed in previous publications. The first criterion is a 
change in sign of the temperature derivative of the resis- 



tivity, dp/dT (see, e.g., Ref. M). A positive (negative) 
sign of the derivative at the lowest achievable tempera- 
tures is empirically associated with a metallic (insulating) 
phase. The second criterion is based on a vanishing acti- 
vation energy, Ea, combined with a vanishing nonlinear- 
ity (threshold voltage, Vc) of current-voltage {I—V) char- 
acteristics when extrapolated from the insulating phase 
|0,||. The activation energy and the threshold volt- 
age are connected via the localization length which is 
temperature-independent. While the derivative method 
deals with the vicinity of the MIT in which the depen- 
dence p{T) is relatively weak, the I — V method is related 
to the insulating phase with exponential p(T). These two 
methods have not been applied simultaneously to the 2D 
MIT. 

In this Letter, we compare these two independent 
criteria for determining the metal-insulator transition 
point in a 2D electron system in silicon metal-oxide- 
semiconductor fleld-effect transistors (MOSFETs). We 
report that in zero magnetic fleld, both methods yield 
the same critical electron density nc(0). Since one of 
the methods is temperature-independent, this equiva- 
lence supports the existence of a true T = MIT in zero 
magnetic field. In contrast, in high parallel magnetic 
fields, where the 2D electrons are fully spin-polarized [||, 
only the I — V method can be used, and it yields the 
critical density ~ 1.5nc(0); the derivative criterion does 
not yield a critical point. This makes uncertain the exis- 
tence of a true metal-insulator transition in a system of 
spin-polarized electrons. 

Measurements were performed on high-mobility silicon 
MOSFETs (^poak = 2.4 x 10* cmVVs at 4.2 K) similar 
to those previously used in Ref. H| . Data were taken by 
a four-terminal dc technique using an electrometer with 
symmetric inputs. 

The resistivity as a function of temperature in zero 
magnetic field is shown in Fig. |l|(a) for several electron 
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FIG. 1. Temperature dependence of the resistivity of a Si MOSFET at different electron densities near the MIT in zero magnetic field 
(a), and in parallel magnetic fields of 4 Tesla (b) and 10 Tesla (c). The electron densities are indicated in units of 10 cm . The inset 
in (c) shows a close-up view of the curve marked by the rectangle. 



densities on both sides of the metal-insulator transition. 
The resistivity of the middle (dashed) curve shows vir- 
tually no temperature dependence over a wide range of 
temperatures; this curve separates those with positive 
and negative dp/dT nearly symmetrically (at tempera- 
tures above 0.2 K) as reported earlier fin]. The existence 
of such a T-independent curve at temperatures down to 
30 mK clearly shows that the logarithmic corrections to 
the resistance (which are expected to be very strong in 
2D systems with resistivity ^ hje^ [p| ) are absent in this 
system in zero magnetic field, see also Refs. 0,p|. Assum- 
ing that the middle curve remains flat down to T = 0, we 
obtain the MIT critical point at n^ = 0.795 x 10" cm"^ 
which corresponds to a resistivity p k, 'ih/e^ , as in other 
experiments on Si MOSFETs 0. We designate the cor- 
responding electron density as the critical density, rici. 

The MIT point has also been determined by studying 
the behavior of nonlinear I — V characteristics on the 
insulating side of the transition. Deep in the insulating 
state {us < rici), a typical low-temperature I — V curve 
is close to a step-like function: the voltage rises abruptly 
at low current and then saturates, as shown in Fig. 0(a). 
The magnitude of the step is 2Vc- (At higher temper- 
atures the curve becomes less sharp, yet the threshold 
voltage, Vc, remains practically unchanged.) Closer to 
the MIT, the I ~ V curves still show a nonlinear step- 
like behavior provided that n^ < rici (see the curve in 
Fig. ||(b) corresponding to n^ = 0.743 x 10^^ cm~^). Ex- 
actly at Us — Uci — 0.795 x 10^^ cm~^, the I — V curve 
is strictly linear (Fig. 0(b)); at Ug > rid, it becomes su- 
perlinear which is characteristic of the metallic state [ [ll| . 
Figure || (closed circles) shows that the square root of the 
threshold voltage is a linear function of electron density 

1/2 

(discussed below). Extrapolation of the Vc {us) depen- 
dence to zero threshold value yields the critical electron 
density ric2- Note that nc2 is equal to rid with high ac- 
curacy. 



Not too deep on the insulating side of the transition, 
and at not too low temperatures, the resistance has an ac- 
tivated form (Tp^l , as shown in the inset of Fig. ||(a) [Q . 
Such a form of temperature dependence was interpreted 
in Ref. ||lj] as a result of the thermal activation of carri- 
ers to the mobility edge, Ec', in this case, the activation 
energy is Ea — Ec — Ep. Since in the insulating regime, 
the I — V curves are strongly nonlinear, we determine 
the resistivity from dV/AI in the linear interval oi I ~V 
curves, i.e., at / ^ 0. Figure S shows Ea as a function of 
the electron density (open circles). In previous detailed 
studies R|, this dependence was found to be linear near 
the mobility edge ||l^. The present data can also be ap- 
proximated by a linear function which yields, within the 
experimental uncertainty, the same critical electron den- 
sity nc2 as the square root of the threshold voltage. The 
reciprocal slope of Ea{ns), D* , can be interpreted (l4| as 
the thermodynamic density of states near the transition 
point. 

The threshold behavior of the I — V curves has been 
explained within the concept of the breakdown in the 
insulating phase (7|jl^. Here we simply outline this con- 
cept. The breakdown occurs when the localized electrons 
at the Fermi level gain enough energy to reach the mobil- 
ity edge in an electric field, Vc/d, over a distance given by 
the localization length, L: eVcL/d = Ec — Ep, where d is 
the distance between the potential probes. The values Ea 
and Vc are related through the localization length which 
is temperature-independent and diverges near the tran- 
sition as L{Ep) oc {Ec — Ep)~^ with exponent s close 
to unity m. This corresponds to a linear dependence 
Vc (us) near the MIT, as seen in Fig. 0. 

We stress that in zero magnetic field, both methods 
— the one based on extrapolation of p{T) to zero tem- 
perature and the other based on the behavior of the 
temperature-independent localization length — give the 
same critical electron density. This adds confidence that 
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FIG. 2. Current- voltage characteristics at a temperature of ~ 
30 iiiK in zero magnetic field. (The aspect ratio is equal to 2.4 for 
this sample.) In case (a), the / — V curve obtained at a higher tem- 
perature (211 niK; dashed line) is also shown for comparison; note 
that the threshold voltage is practically temperature-independent. 
An Arrhenius plot of the resistivity in the insulating phase is dis- 
played in the inset for the following values of B and 71^. T, 
0.741 X 10" cm~2 (circles); 1 T, 0.810 X lO" cm'^ (squares); 
6 T, 0.870 X 10" cm~2 (diamonds). 
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FIG. 3. Activation energy and square root of the threshold voltage 
as a function of electron density in zero magnetic field (circles) and 
in a parallel magnetic field of 4 Tesla (diamonds) . 

the curve with zero derivative dp/dT wiU remain "flat" 
(or at least wiU retain finite resistivity vahie) down to 
zero temperature. 

It is interesting to compare the _B = case with that 
in the presence of a parallel magnetic field. With increas- 
ing parallel field, B, the MIT point nc2, determined from 
the vanishing nonlinearity and activation energy, shifts 
to higher electron densities, saturating above a critical 
field, -Bsat, at a constant value which is approximately 1.5 
times higher than that in zero field (see dots in Fig. 0). 



A similar suppression of the metallic behavior was ob- 
served using a resistance cut-off criterion at the level on 
the order of h/e"^ |l^ . As was shown in Ref. |^ , in the 
metallic phase the saturation of the resistance with par- 
allel field signals the onset of full spin polarization of the 
2D electrons. Hence, one expects that the 2D system is 
spin-polarized at J5 > Bsat, and that the observed phase 
boundary shift is a spin effect. 

Before discussing the temperature dependence of the 
resistivity at B > Bgat, we note that one cannot assume 
that the metallic phase is necessarily strictly character- 
ized by positive dp/dT |18| : one may have a weakly T- 
dependent p{T) with dp/dT < and still have a finite 
resistivity at T = 0. The I—V method yields the electron 
density nc2 at which the exponential divergence of the re- 
sistivity ends, although in principle dp/dT may remain 
negative at this density. 

In Fig. yj(b), we show the temperature dependence of 
the resistivity in a parallel magnetic field high enough to 
cause full spin polarization {B = A Tesla). The mid- 
dle curve corresponds to the critical electron density, 
nc2{B > Bsnt), determined by the method of vanishing 
nonlinearity and activation energy (as shown in Fig. 
by diamonds). In sharp contrast with the B = sit- 
uation, not only are the p(T) curves in the field non- 
symmetric about the middle curve, but all of them have 
negative "insulating-like" derivatives dp/dT < in the 
entire temperature range, although the values of the re- 
sistivity are comparable to those in the B = case. 
Moreover, in a strong parallel magnetic field, there is 
no temperature-independent p{T) curve at any electron 
density: as shown in Fig. ^(c), the curve at the consider- 
ably higher density Us — 1.335 x 10^^ cm~^ compared to 
nc2 — 1.155 X 10^^ cm~^, which could be approximately 
identified as a flat one in the temperature range used, 
changes its slope from weakly- metallic at T > 0.5 K to 
weakly- insulating at lower temperatures. The metallic 
behavior of the resistance as a function of temperature, 
seen at yet higher electron densities in a parallel magnetic 
field (Fig. |l|(c)), is much weaker than in the absence of 
field. We therefore conclude that the derivative method 
does not yield a critical density for the spin-polarized 2D 
system. Its failure leaves uncertain the existence of a true 
metal-insulator transition in a parallel magnetic field. 

Some theories claim that spin-polarized and unpolar- 
ized states are very similar |19[ . The authors of Ref. pc| ] 
considered the temperature-dependent screening of a ran- 
dom potential and predicted metallic {dp/dT > 0) tem- 
perature dependences of the resistivity for both polarized 
and unpolarized states. Therefore, one might expect a 
more or less analogous behavior of p{T) curves around 
the transition in the two cases. However, this is in con- 
tradiction with the experiment: while in zero magnetic 
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The fact that the spin-polarized and spin-unpolarized di- 
lute 2D electron systems behave qualitatively differently 
poses important constraints on the theory. 
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FIG. 4. Dependence of the critical electron density 71c2 (dots) and 
of the thermodynamic density of states at the MIT (squares) on 
parallel magnetic field. The dashed lines are guides to the eye. 

field, "metallic" and "insulating" p{T) curves are approx- 
imately symmetric on both sides of the transition (see 
above), this symmetry completely disappears in a paral- 
lel magnetic field (c/. Figs. |l|(a) and[l|b,c)). An impor- 
tant disagreement between theory |pO| and experiment 
is the dramatic weakening of the metallic temperature 
dependences in the magnetic field, while theoretically, 
the derivative dp/ AT for the spin-polarized state is ex- 
pected to be twice as high for the same electron density 
pO| (for more on this discrepancy, see Ref. |2|]). Thus, 
the properties of the polarized state cannot be deduced 
from those of the unpolarized state in a straightforward 
way Hi. 

Another experimental fact that indicates the effect of 
polarization is the change in the slope of the dependence 
of Ea on Us (Fig. H) which we link, following Ref. p4[ , 
to the thermodynamic density of states at the MIT. In 
zero magnetic field the inverse slope, D*, is close to 
the zero-field density of states, Dt^: D* « 1.2 Dq where 
Dq = 2m/iTh , m = 0.19 TOe, and me is the free elec- 
tron mass. In a parallel magnetic field, D* increases and 
saturates at about 2.7 Dq (see squares in Fig. ^). This 
increase is quite surprising since in the case of full spin 
polarization, one expects that the density of states should 
decrease by a factor of two due to lifting of the spin de- 
generacy |23| ]. So, in many respects the behavior of the 
spin-polarized electron system is peculiar. 

In summary, we have compared two principal ap- 
proaches for determining the critical density for the 
metal-insulator transition in 2D. In zero magnetic field, 
both definitions of the critical point are found to be 
equivalent, strongly supporting the existence of a true 
B — MIT. With increasing parallel magnetic field, the 
I — V criterion gives the critical point which shifts to 
higher electron densities and then saturates, which is 
likely to be a consequence of the spin polarization of 
the 2D electrons. It is accompanied by the disappear- 
ance/weakening of metallic temperature behavior of the 
resistance so that the derivative criterion cannot be used. 
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